HU-TFT-95-71, HD-THEP-95-52, tiep-ph/9512202 



On Bubble Growth and Droplet Decay in Cosmological Phase Transitions 

H. Kurki-Suonio'^'* and M. Laine''''=''l' 
"'Research Institute for Theoretical Physics and ^Department of Physics, 
P.O. Box 9, FIN-OOOU University of Helsinki, Finland 
"Institut fiir Theoretische Physik, Philosophenweg 16, D-69120 Heidelberg, Germany 

(December 1, 1995) 



in 

OS 

o 

Q 



> 

O 
(N 
(N 

in 

OS 



Oh 
I 

Oh 

CD , 
^ ■ 



X 



We study spherically symmetric bubble growth and droplet 
decay in first order cosmological phase transitions, using a nu- 
merical code including both the complete hydrodynamics of 
the problem and a phenomenological model for the micro- 
scopic entropy producing mechanism at the phase transition 
surface. The small-scale effects of finite wall width and sur- 
face tension are thus consistently incorporated. We verify 
the existence of the different hydrodynamical growth modes 
proposed recently and investigate the problem of a decaying 
quark droplet in the QCD phase transition. We find that 
the decaying droplet leaves behind no rarefaction wave, so 
that any baryon number inhomogeneity generated previously 
should survive the decay. 

PACS numbers: 98.80. Cq, 47.75. -l-f, 95.30.Lz 



I. INTRODUCTION 

The cosmological quark-hadron (QCD) phase transi- 
tion took place when the age of the universe was 
t ^ lO^^s. The order of this transition remains undeter- 
mined. If it is of first order, it may take up appreciable 
time. Unless the latent heat is very small compared with 
the surface tension, the probable sequence of events is: 
1 
2 
3 
4 
5 



Supercooling of the quark-gluon plasma. 
Nucleation of bubbles of the hadron gas. 
Rapid growth of hadron bubbles by deflagration. 
Coupling in of particles without strong interactions. 
Collisions of bubbles with shocks from other bubbles, 
and reheating of the plasma to the critical temperature. 

6. Slow growth of the hadron bubbles paced by the ex- 
pansion of the universe, leading to percolation. 

7. Slow shrinking of remaining regions of quark plasma 
and their becoming spherical. 

8. Decoupling of particles without strong interactions. 

9. Evaporation of the quark droplets. 

The electroweak (EW) phase transition |^ took place 
earlier, aXt ^ 10~^^s. It is likely to be of first order, since 
there is baryon asymmetry in the universe However, 
within the Standard Model, the order remains unclear 
for Higgs masses niH ^ 80 GeV. If the transition is of 
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first order, it will likewise proceed through growing bub- 
bles and shrinking droplets, although the details will be 
different. For instance, there may not be a stage of slow 
growth (6-7). 

Clearly, to understand the possible physical conse- 
quences of the cosmological phase transitions, all the 
stages should be thoroughly understood. The nucleation 
period (1-2) from the critical temperature Tc to the nu- 
cleation temperature T„ seems to be rather well under 
control, and we do not touch it here. For the QCD phase 
transition, the present understanding of the other stages 
is the following: The general hydrodynamics of bubble 
growth has been studied extensively Hydrody- 
namics alone cannot determine the bubble solution, but 
due to the strong coupling constant, there is not very 
much more one can say [^Q. The coupling in of elec- 
tromagnetic radiation with the strongly interacting mat- 
ter at about the scale 10^ fm is discussed in |l9|j20f| , the 
collisions with shocks in 0,|8| , and the slow stages at Tc 
in . In our opinion, the possible formation of a shrink- 
ing similarity solution |2^ at stage (7) is not completely 
understood. Assuming that this happens, the decoupling 
of electromagnetic radiation at the scale 10* fm, and the 
final stages until the quark droplet has a radius of about 1 
fm, have been discussed in ||2^,|2^. One should note that 
depending on the values of the latent heat and surface 
tension, the phase transition might also proceed rather 
differently from the standard scenario assumed here . 

In this paper, we complete the hydrodynamic study 
of growing bubbles (3) presented in , and investigate 
in detail the final stages (9) from a radius of about 200 
fm until the time that the quark droplet has completely 
disappeared, thus improving on [2l| . The final stages 
are important since they determine the baryon number 
inhomogeneity resulting from the QCD phase transition. 

For the EW phase transition, there are a number of mi- 
croscopic investigations of bubble growth [p^-pgf . These 
assume that the change of temperature across the phase 
transition surface, 6T/(Tc — Tn), is vanishing, and hence 
are applicable only for very weak deflagrations. A phe- 
nomenological investigation with one free parameter but 
without any limitations on 5T/{Tc — Tn) is in Re- 
cently there has appeared also a microscopic investiga- 
tion with hydrodynamics included [|o|, so that stage (3) 
seems to be reasonably well under control. As to the later 
stages, in ||l^] it was argued that the collisions with the 
shock fronts (5) might have interesting effects, if the bub- 
ble walls are slowed down appreciably. This could happen 
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if the temperature rises close to Tc- Based on simple per- 
turbative estimates, one often assumes that reheating to 
Tc and the stages (6-7) do not exist; however, this need 
not necessarily be the case, see e.g. js^]. In particular, 
lattice simulations suggest that the surface tension 
may be smaller than the perturbative value, leading to 
smaller supercooling and possibly to a reheating to Tc- 
Then the final stages would proceed similarly to the QCD 
case. The results of this paper are applicable to the gen- 
eral characteristics of bubble growth (3) and to the final 
stages (9). 

We study the growth of bubbles and the decay of 
droplets in first-order phase transitions as a hydrody- 
namical problem. We shall use the quark/hadron ter- 
minology, but as stated, our results are qualitative and 
apply to some extent also to the electroweak phase tran- 
sition. A major uncertainty in these studies is the phase 
boundary, the bubble or droplet wall, which remains 
poorly understood, especially for QCD. Hence one usu- 
ally has to resort to some kind of parametrization of its 
properties. We use a cosmic- fluid-order-parameter-field 
model, which allows simultaneous treatment of the phase 
boundary and surrounding hydrodynamics within a sin- 
gle framework p^ ]. We have written a hydrodynami- 
cal code using this model; an earlier version was plane- 
symmetric, but we now have a spherically symmetric ver- 
sion. In the next Section, we describe the model and how 
it differs from the one in [^,^ . In Section III we inves- 



tigate bubble growth with our model, and in Section IV 
droplet decay. 



II. THE MODEL 

We summarize here our cosmic-fluid-order-parameter- 
field model |l8) . We ignore the chemical potential related 
to the small baryon number. The local state of mat- 
ter is then described by three quantities: a local tem- 
perature T, a local flow velocity and a local value 
for the order parameter (f). The order parameter has a 
temperature-dependent effective potential V{(p,T). The 
(meta)stable states of the system are defined by the min- 
ima of the effective potential. 

The equation of state is 



e(0, T) = 3aT^ 
p{c^,T)^aT^-V[<i>,T) 



dV 

V{cp,T)-T—, 



(1) 
(2) 



where a = (7r^/90).g,. The energy-momentum tensor is 

T^"' = a'^0(9> - 5^"'(iaa09"0) + wu^'u'' + g^'^p, (3) 

where w = e + p and the metric convention is ( — 

The energy-momentum conservation equation Tf^^ = is 

split into two parts: 



T^-- (field) = - —CP- = riuW,^^ , 



dV 



dV_ 



(4) 



T.^"^ (fluid) = {wu'^un.t.+p-'' + -g^fb'" - -77M^0,^0'^ (5) 

where ";" is a covariant derivative. Here ?7 is a dissipa- 
tive constant relating entropy production to the gradients 
of (p through T{su'^).fj_ = ri{u^(p^^Y. Eq. (||) contains the 
equation for the (/)-field. 



dV 



(6) 



Eq. (^ (or its equivalent) could in principle be derived 
from field theory at least for the EW case, see e.g. [^ , 
but we do not here attempt to do so. 

In the spherically symmetric case, eq. (^ contains two 
independent equations. For numerical hydrodynamics it 
is best to choose eq. (||) with v = together with the 
contraction of eq. (|^) with u^: 

BV 

{eu^.^+pu^^ - -g^W^^^^ = Viu''4',^?. (7) 

From eqs. (||), (^, (0), the final spherically symmetric 
equations in Minkowski space then become 

- + ^a^(^2g^^) - = VliM + Vdr^), (8) 

dtE + ^dr{r''Ev)+p[da + \dr{r^-iv)\ 
dV 



-rrrliM + vdr^) = vridtcl) + vdr^y, (9) 



1 



dV . 



dtZ + —drir'Zv) + drP + ^d, 

0(p 

= -rijidtcf) + vdr4>)dr(t>, 



(10) 



where E = ej and Z = wj'^v. 

For the effective potential V{4>,T) we use the simple 
parametrization 

y(0, T) = Ij{T' - TS)ct>' - laTct>' + \\c^^ (11) 



where mli] 
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^l + 6a/{Ll,y 
1 

L + 6a/lc 
1 



3aE 



(12) 
(13) 

(14) 
(15) 



Here a is the surface tension, Ic is the correlation length, 
L is the latent heat, and Tc is the critical temperature of 
the transition. We use the functional form in eq. (O) due 
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to its simplicity; our results are qualitative and we expect 
them to remain the same for any V{4i, T) including a first 
order phase transition. 

Let us now compare our model with that used in [ pT|j22| 
(RMPjRM) for studying droplet decay. The main differ- 
ence is that in our model the phase transition surface 
has a microscopic structure and a finite width, whereas 
in ||2l| , p^ it is a discontinuity, across which explicit jump 
conditions are imposed. For this reason RMP have to 
stop the evolution of the decaying droplet at a radius 
of about 1 fm when the width of the interface starts to 
have significance, whereas we are able to follow the decay 
until the very end and beyond. We also believe that in 
our code the velocity determination by the coefficient 
is more natural: we do not need to know anything about 
the structure of the solution beforehand; a fact which, for 
instance, allowed us to find new kinds of solutions in [p^ . 
In 11^,^ , the velocity is fixed by giving the ratio of the 
hydrodynamical and thermal fluxes at the deflagration 
front. 

On the other hand, the code in |2|,|2^ contains fea- 
tures we do not have. First, RMP have an expanding 
background metric. We do not expect the extremely slow 
expansion to have any effect at the short time scales we 
are investigating. Second, in |^ the electromagnetic ra- 
diation is treated properly, allowing RM to investigate 
the decoupling of radiation from the strongly interacting 
matter at a radius of about 10'* fm, which is important for 
concentrating baryon number inside the quark droplet. 
Here we study only the final stages when the bubble has 
a radius smaller then 200 fm, so that the decoupling has 
already taken place. Hence our number of degrees of 
freedom is that of the strongly interacting particles. We 
take the physical result of the decoupling period, namely 
a strongly increased baryon number density within a ra- 
dius of 10^ — 10'' fm around the droplet |Q, as an initial 
condition for our analysis of droplet decay. 



III. BUBBLE GROWTH MECHANISMS 



does not forbid any of the solutions, and the parameter 
r] is allowed to vary freely. At least for the QCD case, 
there is ample parameter space for all three kinds of so- 
lutions m. 

We set out to verify that all three classes can be re- 
alized. We choose parameter values which allow the 
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0.1T3, I, = 6T- 



three classes: L 
a = 34.257r^/90. The nucleation temperature corre- 
sponding to these is Tjnit = 0.86Tc Q. We start with 
initial data where there is a small "recently nucleated" 
bubble surrounded by homogeneous fluid at rest. The 
bubble will start to grow and we follow it until the con- 
figuration becomes self-similar, see Fig. ^ 
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FIG. 1. The velocity profile approaches the similarity so- 
lution as the bubble grows. This run is for rj = l.OTc. The 
time interval between the dashed curves is At — 15S.75T^^ , 
and between the solid curves At = 675r<r\ The horizontal 
axis is ^ = r/t, distance scaled with time. 



In |l6| we discussed the different hydrodynamical 
growth mechanisms of phase transition bubbles. We 
found three different classes: (1) weak deflagrations mov- 
ing at a subsonic speed and preceded by a shock front and 
a compression wave; (2) Jouguet deflagrations moving at 
a supersonic speed, preceded by a shock front and a com- 
pression wave, and followed by a rarefaction wave; and 
(3) weak detonations moving at a supersonic speed, and 
followed by a rarefaction wave. 

Note that for small supercooling, the condition of non- 
negative entropy production allows only weak deflagra- 
tions |^,||,^. Moreover, adding information about the 
microscopic entropy producing mechanism, the solution 
becomes fixed; for instance, for the EW case mildly rela- 
tivistic weak deflagrations seem probable Q . The above 
three classes refer to a case when the entropy condition 



The velocity of the bubble wall will depend on the dissi- 
pative constant rj. We covered the range from 77 = O.OlTc 
to 7y = lOTc, and found all three classes: (1) weak defla- 
grations for rj > 0.13rc; (2) Jouguet deflagrations for 
0.12rc ^ '7 ^ 0.13rc; and (3) weak detonations for 
rj < 0.12Tc. In Fig. || we show a sequence of final ve- 
locity profiles for a set of runs, and in Fig. ^ we show the 
wall and shock velocities as a function of rj. The shock 
fronting the compression wave in the deflagration bubble 
solutions becomes exponentially weak for slowly growing 
bubbles, and we did not resolve the shock for 77 > O.lTTc. 
Fig. 1^ shows an example from each solution class. 

We see (Fig. ||) that as 77 is decreased the wall velocity 
grows smoothly until it is close to the speed of sound. 
Then the wall velocity appears to have some difficulty 
breaking the sound barrier. (In a Jouguet deflagration 
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FIG. 2. A sequence of velocity profiles for bubbles with 
different values of (given in units of Tc next to each profile) . 
These profiles are from our hydrodynamical runs. Compare 
with Fig. 3 of |l6j , where the profiles are solutions of eqs. ( |2c| ) 
and (|2l[) (with a bag equation of state). 



the wall is subsonic with respect to the fluid just ahead 
of the wall, and it is moving at the speed of sound with 
respect to the fluid just behind, although the velocity 
with respect to the bubble center, "the wall velocity", 
already exceeds the speed of sound). The increase in 
wall velocity slows down, until at 77 ^ 0.12Tc the solution 
shifts from a deflagration to a detonation solution and 
the wall velocity jumps abruptly. Thus we do not have 
solutions for wall velocities between v — 0.593 and v — 
0.635. 

In 1^ we described similarity solutions for all wall 
velocities. The reason that we do not here find solutions 
in the above range 0.593 . . . 0.635, is that there are pairs 
of solutions with different wall velocities corresponding 
to the same value of 77. The slower one is a Jouguet 
deflagration, the faster one a weak detonation. Only one 
of these two is realized for a given initial configuration. 
Thus when the solution shifts from a defiagration to a 
detonation, there is a jump in wall velocity. In fact, in 
the runs with the fastest defiagrations, the bubble first 
went into a detonation configuration before settling into 
a deflagration, see Fig. ||. 



IV. DROPLET DECAY 

Kajantie and Kurki-Suonio (KK) have discussed the 
shrinking quark droplet. There is an outward fluid flow 
from the surface of the droplet. When the droplet van- 



FIG. 3. The bubble wall velocity and the shock velocity as 
a function of rj. 



ishes the source of this flow disappears. This led KK to 
conclude that the evaporation site sends out a rarefac- 
tion wave by which the outward flow is stopped. KK 
did not consider the effect of surface tension on droplet 
evaporation. 

RezzoUa, Miller, and Pantano |^ have improved upon 
the work of KK. They conclude that the ffow pattern 
around a shrinking droplet settles into a similarity solu- 
tion. If the similarity solution would hold till the end, the 
droplet would leave behind a homogeneous ffuid at rest, 
since all structure around the droplet would have shrunk 
to zero size at the moment the droplet vanishes. But 
as the droplet becomes small, the surface tension causes 
the droplet wall to accelerate, and the droplet actually 
vanishes before the surrounding outward flow pattern. 
Apparently RMP have not followed the evolution past 
the droplet disappearance with their code. This we have 
done with ours. 

We start with the similarity solution as an initial con- 
dition. Thus we need to construct this solution in our 
model. Let us first discuss the regions away from the 
phase transition surface, where the order parameter lies 
at the minima (/> = and — 0min(r) in the two phases 
{dV/d(t) = 0). We assume the distance scale is large 
enough so that 9^0 terms can be ignored in T^". Then 
we have 



3aT^ 



p,,(T) = aT4-y(0,,i„(r),r). 



T 



dV 
df 



(16) 
(17) 

(0„,in(T),r), (18) 
(19) 
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FIG. 4. The energy density, velocity, and order parameter 
profiles for the three different solution classes. The first (a) 
is a weak deflagration, the second (b) a Jouguet deflagration, 
and the third (c) a weak detonation. 
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FIG. 5. A bubble (77 = 0.1207Tc) which first goes to a 
detonation configuration before settling to a deflagration. The 
three profiles are from times tx = ISSOT^'S £2 = 2700T7\ 
and tz = 10800r<rV 



The similarity solution depends on the coordinates r 
and t only through the combination ^ = r/t. For a 
shrinking solution, we must choose t = at the mo- 
ment of disappearance, so t < and thus ^ < 0. The 
similarity solution has homogeneous fluid at rest inside 
the droplet, v = 0, e = eq{Tq). The droplet surface at 
£, = ^dofl is a deflagration front, the phase boundary mov- 
ing in at a speed |^dcfl|- The usual deflagration conditions 
give us then the fluid state just outside the droplet wall, 
with fluid flowing outwards. The profile e(^), of this 
compression wave is solved from Mm 



1 



iv-o'-a-<r 

dp 



dv 
dT 



2-(l-<)(l-«2), (20) 



}^ dp ,rj.. dT ^ C-v 1 dv 

e{T)+p{T)dT^ ' di l-vil-v^di' ' 

The velocity is positive (outwards) and becomes zero at 
the sonic point ^ = — c^, where we have a weak discon- 
tinuity. For 1^1 > Cs the fluid is at rest with constant 
energy density. (Note that there is no shock front, un- 
like in the case of an expanding deflagration bubble). In 
Fig. H we show a similarity solution for the case Tq = 0.99, 
Cdcfl = -0.05. 

For the phase transition surface we use a finite- 
width wall obtained with our stationary plane-symmetric 
code This gives the initial data for our dynamical 

code. 

We present results from a run with parameter values 
0.5T3, = ITj-i, and a = 407rV90, using 
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FIG. 6. A similarity solution for a shrinking droplet. The 
two values in the brackets refer to the value of the correspond- 
ing quantity at "0" and "1" of the y-axis. 



an initial condition with Tq = 0.99Tc. (This is a fairly 
arbitrary choice, as there is large uncertainty about the 
values of these parameters, and we are thus looking for 
qualitative features only). The initial data is the simi- 
larity solution for the case Wdofi — —0.05, and we set the 
initial radius of the droplet to be r = 200T-^ (for QCD 
is just 1 fm, and for the EW case 10~^ fm). Now the 
actual wall velocity will depend on the dissipation param- 
eter 77. The value of 77 corresponding to a given Vdcfi can 
be found by trying out different values. Because of the 
fairly small initial size, there are already some small-scale 
effects and the solution should be deviating somewhat 
from the similarity solution. Therefore using the similar- 
ity solution as initial data will cause some disturbance to 
spread out from the wall to alter the profile. We min- 
imized this initial disturbance by choosing rj = 1.97Tc, 
giving an actual initial wall velocity Wdofl — —0.047. 

Figs. |7|-p0| show results from this run. At first we stay 
fairly close to the similarity solution but as the droplet 
gets smaller the wall is accelerated inwards, and the tem- 
perature and energy density rise inside the droplet. As 
the droplet disappears the outward flow leads to an un- 
derdense region with lower temperature in the center. 
The fluid flow then turns around filling this underdense 
region. A pulse of inward flow is thus sent outwards, and 
decays rapidly. The end result is homogeneous fluid at 
rest. Thus we see no rarefaction wave propagating to in- 
finity after the decay. This is in contrast to expectations 

in 

In Fig7l9 we show the compression factor p (density of 
noninteracting test particles moving with the fluid). In 



the end the thermodynamical properties {T,v,(j>) of the 
fluid become homogeneous. There is a small depression 
in the final p at the droplet evaporation site. This is 
due to the droplet surface energy which was converted 
into thermal energy. Because of this extra energy some 
decompression was needed to even out the energy density. 

Note that due to the difficulty of getting baryon num- 
ber through the wall, the baryon number tends to concen- 
trate on the inside of the droplet surface, and is thus not 
moving with the fluid [p],p3|-p7 22 1. Thus the compres- 
sion factor p does not represent baryon number density. 
However, it can be viewed as representing a background 
motion of the baryon number, upon which the baryon- 
number-concentrating effects are superimposed. Thus we 
see that the hydrodynamical flow itself does not lead to 
any inhomogeneity in the baryon number, except the tiny 
depression at the droplet evaporation site. 

RM 1 22 1 studied the effect of radiative energy transfer 
out of the quark droplet on the baryon number distri- 
bution. We are studying scales where this transfer has 
ceased to operate as the droplet has become transpar- 
ent to radiation. Thus we can combine our work with 
RM by considering their result at the relevant time as 
an initial condition for our runs. Their compression fac- 
tor prm can be viewed as being proportional to baryon 
number density in an idealization where other baryon 
number concentrating effects are ignored. Let us do this 
idealization in the discussion below. RM found that the 
shrinking droplet leaves behind a region with a high over- 
density in baryon number, whose radius is of the order of 
the mean free path A ^ lO^fm of the electromagnetic ra- 
diation. Inside the quark droplet there is an even higher 
baryon number density. As we are considering the region 
r ^ A, we can take our compression factor p to denote 
the baryon number density normalized to the value it has 
in the end at about r ^ lO'^fm, where it is fairly homoge- 
nous. Thus our Fig. |^(a) shows the baryon number den- 
sity contrast between the quark and hadron phases (the 
numerical value is much smaller in our case than in |2^ ] 
because of different parameter values used). There are 
two contributions to this density contrast: (1) the de- 
compression due to the change in equation of state as 
the fluid passes through the phase boundary, and (2) the 
extra overdensity inside the droplet due to surface ten- 
sion, which becomes significant in the end. There will 
be no remnant of this higher density inside the droplet 
after the droplet has disappeared (in reality, there prob- 
ably will be a remnant due to effects we are ignoring 
here |^5|). Contribution (1) is eliminated as the fluid is 
converted to the hadron phase. Contribution (2) is elim- 
inated as the thermodynamical variables are evened out 
by the outward pulse (see Fig. ||). Thus the /Orm outside 
the droplet represents the final baryon number distribu- 
tion, and it is not disturbed by any rarefaction wave. 
For r <C A the baryon number density is homogeneous 
except for the central small depression caused by the re- 
leased droplet surface energy (see Fig. ^(b)). The baryon 
number overdensity within r ^ lO^fm is later eliminated 
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mainly by neutron diffusion psf . 
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V. CONCLUSIONS 

We have studied the growth of bubbles and decay of 
droplets in cosmological phase transitions using a spheri- 
cally symmetric version of our hydrodynamic code based 
on a cosmic- fluid-order-parameter model 1 18| of the tran- 
sition. 

We demonstrate how an initial small newly nucleated 
bubble begins to grow and evolves to a similarity solution. 
The wall velocity depends on the dissipative constant 77. 
We find all three classes of hydrodynamic similarity so- 
lutions: weak deflagrations, Jouguet deflagrations, and 
weak detonations. In the first class the wall velocity is 
subsonic, in the other two classes supersonic with respect 
to the origin. 

Not all wall velocities are realized as the value of 77 is 
varied. For a certain range in rj, there would be both 
a Jouguet deflagration and a weak detonation solution 
with the same value of 77, but a different wall velocity. 
As the solution shifts from the former class to the latter, 
there is thus a jump in wall velocity. 

Note that all classes of solutions are not possible for all 
values of the other parameters [||j|j2^. Typically weak 
deflagrations are always possible, but detonations and 
Jouguet deflagrations may not be reached; or the solution 
may shift directly from weak deflagrations to weak det- 
onations as 77 is decreased, skipping the class of Jouguet 
deflagrations [ p3| . 

We followed the evaporation of a quark droplet start- 
ing from a similarity solution. The droplet deviates from 
the similarity solution in the end when small-scale effects 
from the surface tension and finite wall thickness become 
significant. The nature of the fluid flow after the droplet 
evaporation differs from what has been previously sug- 
gested |^j2l|,^. The deviation from similarity flow in 
the end leads to a temporary energy overdensity at the 
evaporation site, surrounded by a shell of underdensity. 
This inhomogeneity is eliminated by an outward moving 
pulse which decays rapidly. This effect does not extend 
beyond the small scales where the droplet deviated from 
the similarity solution. 

Our model does not contain the effects which concen- 
trate baryon number: the suppression of baryon flow 
through the wall and the radiative entropy transfer out of 
the droplet p|p^- p6| , ^ . The temporary hydrodynamic 
compression at the final stage of droplet decay does not 
affect the baryon/entropy ratio and thus has no lasting 
effect on the baryon density. There is a small decom- 
pression at the very center to accommodate the droplet 
surface energy which has been converted to thermal en- 
ergy. The absence of a global rarefaction wave indicates 
that any previously generated baryon number inhomo- 
geneity at the scale 10'' fm is not diluted away by the 
hydrodynamics related to droplet decay. 
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FIG. 7. Results from the run. We show profiles for four 
quantities: proper energy density (a), local temperature (b), 
order parameter (c), and flow velocity (d). The solid curves 
are from successive time-slices as the droplet is shrinking (wall 
moving to the left). The dashed curves are from time-slices af- 
ter the droplet has disappeared (pattern moving to the right) . 
We use a logarithmic axis for the radial coordinate to focus 
on the droplet evaporation site. The finite thickness of the 
wall becomes thus apparent when the wall has reached small 
values of r. 



FIG. 8. Results from the same run as Fig. |^. We show 
later time slices (after droplet evaporation) of proper energy 
density (a) and flow velocity (b) with an expanded scale. A 
pulse of inward flow is moving outwards and decaying rapidly, 
evening out the energy density (and other thermodynamic 
quantities). 
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FIG. 9. The compression factor p. The upper figure (a) 
is like Fig. |^ The lower figure (b) is a close-up near the 
droplet evaporation site and shows time slices after droplet 
evaporation. The solid line is the final p-profile, showing the 
decompression due to droplet surface energy. Note the very 
small spatial extension of the decompressed region. 
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FIG. 10. The flow lines of particles moving with the fluid. 
The lower figure (b) is a magnification of the final stages. This 
figure contains the same information as Fig. ^ but in a dif- 
ferent form. At the time of the droplet decay, there is a local 
rarefaction and subsequent compression near the droplet, but 
no global rarefaction wave is left behind. Due to the difficulty 
of getting baryon number across the phase transition surface 
(the thick line), a part of the baryon number would not follow 
the flow lines but would instead remain in the quark phase. 
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